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Abstract. This paper is in the form of an essay. It defines fractal tops and 
code space structures associated with set-attractors of hyperbolic iterated func- 
tion systems (IFSs). The fractal top of an IFS is associated with a certain shift 
invariant subspace of code space, whence the entropy of the IFS, and of its 
set-attractor, may be defined. Given any ordered pair of hyperbolic IFSs, each 
with the same number of maps, there is a natural transformation, constructed 
with the aid of fractal tops, whose domain is the attractor Aiof the first IFS 
and whose range is contained in the attractor A2 of the second IFS. This trans- 
formation is continuous when the code space structure of the IFS is "contained 
in" the code space structure of the second IFS, and is a homeomorphism be- 
tween A\ and A2 when the code space structures are the same. Conversely, 
if two IFS are homcomorphic then they possess the same code space struc- 
ture. Hence we obtain that two IFS attractors are homeomorphic then they 
have the same entropy. Several examples of fractal transformations and fractal 
homeomphisms are given. 



1. Introduction 

In this essay we introduce fractal transformations. The main examples are fasci- 
nating mappings between diverse subsets of R 2 ; they can be readily illustrated by 
using the chaos game. Fractal transformations can be quicky grasped because they 
rely on basic notions in topology, probability, dynamical systems, and geometry. 
They may be applied to computer graphics to produce digital content with new 
look-and-feel [I]; they may also be relevant to image compression and biological 
modelling. 

2. Hyperbolic IFS 

Definition 1. Let (X, dx) be a complete metric space. Let {/1, f 2 , /n} be a finite 
sequence of strictly contractive transformations, f n : X — > X. for n = 1,2, N. 
Then 

T:={X;f 1 J 2l ...J N } 
is called a hyperbolic iterated function system or hyperbolic IFS. 



Date: March 13 th 2007. 

1991 Mathematics Subject Classification. [2000]Primary 05C38, 15A15; Secondary 05A15, 
15A18. 

Key words and phrases. Iterated function systems, fractal geometry, dynamical systems, in- 
formation theory. 

The author thanks Louisa Barnsley for help with the illustrations. 

This paper is in final form and no version of it will be submitted for publication elsewhere. 

1 



2 



MICHAEL F. BARNSLEY 



A transformation f n : X — * X is strictly contractive iff there exists a number 
l n € [0,1) such that d(f n (x),f n (y)) < l n d(x,y) for all x,y £ X. The number i„ is 
called a contractivity factor for /„ and the number 

I = max{Zi,Z2, In} 

is called a contractivity factor for T . 

Let f2 denote the set of all infinite sequences of symbols {(Jk}kLi belonging to 
the alphabet {1, ...,7V}. We write a = aia 2 a 3 ... £ £1 to denote a typical element of 
fi, and we write u>k to denote the k th element of u> S f2. Then (17, dn) is a compact 
metric space, where the metric da is defined by dn(cr, u>) — when a = lo and 
cfo(<7, w) = 2~ fc when 7c is the least index for which a% ^ u>k- We call fl the code 
space associated with the IFS J- . 

Let a € f2 and a; € X. Then, using the contractivity of .T 7 , it is straightfoward to 
prove that 

<?V(cr) := lim f ai o j~ & 2 a ...f ak {x) 

k — >cjo 

exists, uniformly for x in any fixed compact subset of X, and depends continuously 
on a. See for example [T], Theorem 3. Let 

Af = {</VO) -cr en}. 

Then Ajr C X is called the attractor of T . The continuous function 

is called the address function of T. We call := {<r e il : (f)^(cr) = x} the 

set of addresses of the point x £ A jr. 

Clearly Ayr is compact, nonempty, and has the property 

A r = A(At) U J 2 (At) U ... U f N (A F ). 

Indeed, if we define H(X) to be the set of nonempty compact subsets of X, and we 
define T : H(X) -> H(X) by 

(2.1) F{S) = fi(S) U fc(S) U ... U f N (S), 

for all S 6 H(X), then At can be characterized as the unique fixed point of T , see 
[TT] . section 3.2, and [22] . 

IFSs may be used to represent diverse subsets of M 2 . For example, let A, B, and 
C, denote three noncollinear points in R 2 . Let a denote a point on the line segment 
AB, let b denote a point on the line segment BC and let c denote a point on the 
line segment CA, such that {a, b, c} fl {A, _B, C} = 0, see panel (i) of Figure [TJ Let 
hi : M 2 — > R 2 denote the unique affine transformation such that 

hi(ABC) = aBc, 

by which we mean that h\ maps A to a, B to 75, and C to c. Using the same 
notation, let affine transformations 7i2, ^3, and /14 be uniquely defined by 

h 2 {ABC) = abC, h 3 (ABC) = Abe, and h±(ABC) = abc. 

Let Tafin = {R 2 ; hi, h 2 , h%, /14} where a = |£?c|/|AB|,/? = \Ca\/\BC\ and 7 = 
j46|/|C/1|. The attractor of T a ,p,-i is the filled triangle with vertices at A, B, 
and C. The attractor of the IFS {R 2 ; hi, h 2 , /13} is an affine Sierpinski triangle, as 
illustrated in (iii) in Figure [TJ 
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(i) (ii) (iii) 

Figure 1. (i) The points used to define the affine transformations 
h n : R 2 — > R 2 for n = 1,2,3,4; (ii) sketch of the attractor of 
the IFS {R 2 ; hi, hi, /13, hi}; (iii) sketch of the attractor of the IFS 
{R 2 ; hiMM}- Here a = 0.65, (3 = 0.3, and 7 = 0.4. 
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For reference we note that when A = (0,0), B = (0,1), and C — (0.5,1), the 
transformations of the IFS T a fi^ are given by 

h„(x, y) = (a n x + b n y + c„, d n x + e n y + l n ) 

with the parameters specified in Table [TJ We will write ▲ to denote the filled 
triangle ABC. 

3. Chaos game 

When the underlying space is the euclidean plane, one way to sketch the attactor 
Ayr of an IFS T is to plot the set of points 

U = {f ai o / CT2 o ...f aK (x) : a k e {1, 2, N}, k = 1, 2, ...if}, 

for some 1 6 X and some integer if. The Hausdorff distance between Ajr and Ajf 
is bounded above by C ■ l K where the constant C depends only on T and x. 

A more efficient method is by means of a type of Markov Chain Monte Carlo 
algorithm which we refer to as the chaos game. Starting from any point (xq, yo) € 
R 2 , a sequence of a million or more points {(xk, Uk)}^ = o is computed recursively; 
at the k th iteration one of the functions of T is chosen at random, independently of 
all other choices, and applied to (xk-i, J/fe-i) to produce (xk,yk) which is plotted 
when k > 100. The result will be usually a sketch of the attractor of the IFS, 
accurate to within viewing resolution. 



4 



MICHAEL F. BARNSLEY 



The reason that the chaos game yields, almost always, a "picture" of the attractor 
of an IFS depends on Birkhoff's ergodic theorem, see for example [9]. The scholarly 
history of the chaos game is discussed in [T2] and [3T], and appears to begin in 1935 
with the work of Onicescu and Mihok, [17]. Mandelbrot used a version of it to 
help compute pictures of certain Julia sets, [15] pp. 196-199; it was introduced to 
IFS theory and developed by the author and coworkers, see for example [I], [2], 
[6], and [8], where the relevant theorems and much discussion can be found. Its 
applications to fractal geometry were popularized initially by the author and others, 
see for example Q], [7J, [T5], and [20] , 

The sketches in panels (ii) and (iii) of Figure [1] were computed using the chaos 
game. At each iteration the function h n was selected with probability proportional 
to the area of the triangle h n (ABC), for n = 1,2, 3, 4. 

In section [JJ we show how the chaos game may be modified to calculate examples 
of the fractal transformations that are the subject of this article. Hopefully you 
will be inspired to try this new application of the chaos game. 

4. The tops function 
We order the elements of il according to 

a < lo iff o-fe > uj k 

where k is the least index for which o~k ^ u>k- This is a linear ordering, sometimes 
called the lexicographic ordering. 

Notice that all elements of f2 are less than or equal to 1 = 11111... and greater 

than or equal to N = NNNNN Also, any pair of distinct elements of £1 is 

such that one member of the pair is strictly greater than the other. In particular, 
the set of addresses of a point x G Ayr is both closed and bounded above by 1. It 
follows that <i>yr X (\x\) possesses a unique largest element. We denote this element 
by Tjr[x). 

Definition 2. Let T be a hyperbolic IFS with attractor Ayr and address function 
4>jr : Ci _> Ayr. Let 

Tp{x) = max{er G D, : <f>r{o-) = x} for all x € Ajr . 

Then 

■= {tf{x) : X £ Ayr} 

is called the tops code space and 

Tjr : Ayr — > ttyr 

is called the tops function, for the IFS T . 

Notice that the tops function Tyr is one-to-one. It provides a right-hand inverse 
to the address function, according to 

(j>jrOTjr = i A:F 

where iA? denotes the identity function on Ayr and o denotes composition of func- 
tions. Let 

<$>yr : n r -> Ayr 

denote the restriction of 4>yr to flyr, defined by &yr(a) = <fiyr(a) for all a G ilyr. 
Then <&yr is the inverse of Tyr, namely 
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We note that although <&jf is one-to-one, onto, and continuous, Tyr may not be 
continuous. Let Cljr denote the closure of flyr, treated as a subset of the metric 
space (O, (in)- Let 

$jr : Of -> Ayr 

denote the restriction of 4>yr to SI jr. Then &jr is continuous and onto. Notice that 
the ranges of <i>jc- and &jr are both equal to Ajr because A? is closed. 

5. Fractal transformations 

Let Q denote a hyperbolic IFS that consists of N functions. Then 4>g o tjt : 
Ayr — » Ag is a mapping from the attractor of T into the attractor of Q. We refer 
to 4>g o tjt as a fractal transformation. 

In order to illustrate transformations between subsets of K 2 we use pictures. We 
define a picture function to be a function of the form 

where £ is a color space. A picture function assigns a unique color to each point 
in its domain Dtp. For example we may have €. ={0, 1, ...255} 3 and each point of £ 
may specify the red, green, and blue components of a color. A picture in the non- 
mathematical sense may be thought of as a physical representation of the graph of 
a picture function. 

If T : cl 2 -> then £2 = *}3 o T denotes a picture whose domain is Dq . 
We can obtain insights into the nature of T by comparing the picture functions 
and o T, where represents a given picture which may be varied. We will use 
this method to illustrate fractal transformations. 

Let 

T = {C; fi{z) = sz - 1, h{z) =az + l, for all z G C}, 
Q = {C; g{z) — s±z — 1, g{z) = s±z + 1, for all z € C}, 
H = {C; h(z) = s 2 z - I, h(z) = s 2 z + 1, for all z E C}, 
where C, denotes the complex plane, s = 0.5(1 + i), Si = 0.44(1 + i), and s 2 — 
0.535(1 + i). We denote the attractors of these IFSs by Ajr 7 Ag, and An- In the 
top row of Figure [2] we illustrate, from left to right, the three picture functions, 
*Pe : Ag — » C, : Ajr — > £, and : C These pictures were obtained by 

masking a single original digital picture, whose domain we took to be {z — x + iy G 
C : —3.5 < x < 3.5, —3.5 < y < 3.5}, by the complement of each of the sets Ag, 
Ajr, and An- 

The attractor Ayr is a so-called twin-dragon fractal. It is an example of a just- 
touching attractor: that is, f\{Ajr) n f 2 (Ayr) is non-empty and equals f\{dAjr) n 
f 2 (dAjr) where dAjr denotes the boundary of Ayr. This contrasts with Ag which 
is totally disconnected, perfect, and in fact homeomorphic to the classical cantor 
set. This also contrasts with Au which is such that there exists a disk in R 2 , of 
non-zero radius, which is contained in hi(A-n) H h 2 (A-)-c). 

The bottom row of Figure[2]illustrates the pictures, from left to right, tygocjjgoTyr, 
tyyrocfyyr o Tyr , and < tyn 0< t } TL ° T T ■ They were computed by a variant of the chaos game 
as explained in section [7] The domain of each of these pictures is Ayr. We notice 
that tyyrotfiyroTyr = yr , which is true regardless of the choice of IFS J- since <\>yroryr 
is the identity on Ayr. We notice that both tygocfygoTyr and ^n 0( t ) n 0T J r have features 
in common with the underlying digital picture; for example, tygofig ° Tyr displays 
something like the texture of the hat, near the middle of the bottom-left image. 
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Figure 2 . Examples of fractal transformations are illustrated us- 
ing the three picture functions tyg, and tyn shown in the 
top row, from left to right. The domains of these functions are 
the attractors Ag, Ajr 1 An of three IFSs J 7 , Q, and H, defined 
in the text. The bottom row illustrates the pictures tygocfig o tjf, 

tyjrOtfijr O Tjr, and ^H°4>H ° TjF- 

The bottom right image shows parts of the hat, repeated several times, and some 
clearly delineated small twin-dragon tiles. 

We are led to consider the following questions. Under what conditions on general 
IFSs T and Q is the fractal transformation (figorjr continuous? When does it provide 
a homeomorphism between Ayr and Ag? 

Definition 3. The address structure of T is defined to be the set of sets 



The address structure of an IFS is a certain partition of fljr. Let Cg denote the 
address structure of Q. Let us write C? -< Cg to mean that for each S £ there 
is T G Cg such that S C T. Notice that if Cjr = Cg then fijr = fig. Some examples 
of address structures are given in section [6l 

Theorem 1. Let T and Q be two hyperbolic IFSs such that Cjr ~< Cg. Then the 
fractal transformation <j)g o tjt : Ajr — * Ag is continuous. If Cjr = Cg then (fig o tjt 
is a homeomorphism. 

The proof relies on a standard result in topology, Lemma 2 below, which we 
present in the context of metric spaces. 

Lemma 1. (cf. [16] . bottom of p.194-.) Let F : X — > Y be a continuous mapping 
from a compact metric space X onto a metric space Y . Then S C Y is open if and 
only if F~ 1 (S) C X is open. 



Cjr = {(fi^ 1 {{x})nn^:xeA r }. 
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Proof. If S C Y is open then F^ 1 (S) C X is open because F : X — > Y is continuous. 
Suppose that F _1 (5) is open. Then X\F~ 1 (S) is closed. But a closed subset 
of a compact metric space is compact. The continuity of F now implies that 
F(X\F~ 1 (S)) is compact and hence closed. But F(X\F~ 1 {S)) = Y\S. Hence 
S is open. □ 

Lemma 2. (cf. |16j . Proposition 7.4 on p. 195.) Let F : X — > Y be a continuous 
mapping from a compact metric space X onto a metric space Y . Let H : Y — ► Z 
where Z is a metric space. Let H o F : X — > Z be continuous. Then H : Y — > Z is 
continuous. 

Proof. Let O C Z be open. Then (fl"oF) -1 (0) = F -1 (fl" -1 (0)) is open. But then 
by Lemma 1 i/ _1 (0) is open. Hence H : Y — > Z is continuous. □ 

Proof of Theorem 2. In Lemma 2 we set X = Q,p, Y = Ap, and Z = Ag. We 
choose F : X — > Y" to be $jr : f2jr — > Ajr. Then F : X — * y is a continuous 
mapping from a compact metric space X onto a metric space y. We also choose 
H : Y -> Z to be 

H = 4>g o Tyr : Ayr — y Ag . 

Now look at the function 

G := H o F = <j)g o rjr o $jr : fi^r — > Ag. 

If er £ f2^r, then both cr and (ry o < E > ^-)(ct) belong to the same set in the Cp. Since 
Cjr -< Cg it follows that both a and (ry o $jf)(ct) belong to the same set in Cg. It 
follows that 

{4>g °Tjr o $jt)(ct) = <fig(cr) for all cr £ f2jF. 

But 0g : f2 — > is continuous. Hence G is continuous. 

We have shown that the conditions in Lemma 2 hold. It follows that H = <f>goTf 
is continuous. 

When Cjr = Cg it is readily verified that <f)g o tjp : tLf — » Ag is one-to-one and 
onto and that its inverse is fir o rg. Also = Cg implies Cg -< CV- and so, by the 
first part of the theorem, c/>jf o rg is continuous. Hence 4>g o : Ajr — > Ag is a 
homeomorphism. □ 



6. Examples of address structures 

6.1. Backwards orbits. Let x £ Ap. Let a £ ft be such that 4>r(cr) = x. 
Assume that the /„'s are one-to-one. Then define xq — x and xt — f^Hx^-x) for 
k = 1,2,3, .... Notice that Xk = (j)^(crkO'k+iO'k+2---)- We call {x/c}^ a backwards 
orbit of x (under the IFS J 7 ). 

The set of all addresses of x can be calculated by following all possible back- 
wards orbits of x. Define a sequence of points {xk}kLo m A^- and an address 
a = a{a2<Jz... £ ^l, as follows. Let xq — x. For each k = 1, 2, 3, ... first choose 

a k £ {n £ {1, 2, N} : x k ^ £ /„(Ajr)} 

and then define 

Xk = f~^(x k -i)- 

Then o~ £ i/ijp 1 ({a;}) and all a £ </ , j(r 1 ({2 ; }) can be obtained in this manner. 
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Figure 3. Before, on the left, and after, on the right, a fractal 
homcomorphism. See text. 

6.2. Some notation. We use the notation (PQ) = PQ\{P, Q} to denote the 
straight line segment which connects the two points P and Q in R 2 , without its 
endpoints. We write ft' to denote the set of all finite length strings of symbols 
from the alphabet {1, 2, N}, including the empty string "0" . We write \<r\ 
to denote the length of a G f2'. We define uia — u)\Ll}-2,...lo\ u \0\02---0\ ( ,\ for all 
uj,a G ft'. Similarly we define u>a — u)iu>2—^\u:\< 7 i cr 2--- for all w G ft', a G ft. We 
write S : ft' — ► ft' to denote the shift operator defined by S(<r) — a%az.--cr\ ( j\ when 
\a\ > 1 and S("0") ="0". We write / CT = f ai o / ff2 o ... o for all cr G ft' with 
|ct| > 1 and denotes the identity function. 

6.3. Example 1. An interesting example of address structures is provided by the 
IFS T = J- a .j3 n = {R 2 ; hi, h>2, /14}, introduced at the end of section^ Here we 
prove that 

for all a, /3, 7, a, (3, 7 € (0, 1) by calculating the address structure C?^ p . 

In this case there is only one backwards orbit {xk}^ of each x G Ap. This 
is because of the form of T , and because the h n 's are affine and so preserve ratios 
of distances between points which lie on any given straight line: for example if 
x G ab = /12(A) R /14(A) then /i 2 ~ 1 (x) = h^ 1 (x). Indeed, the mapping T : A — > A 
defined as in equations 18.11 is continuous and the backwards orbit of x is the same 
as the orbit of x under T treated as a dynamical system. 

Any point xk G {xk}^ on the backwards orbit of x, such that more than one 
map h~ x may be applied, is such that xk belongs to the set 

U ftj(A) n h.j{k) = (ab) U (be) U (ca) U {a,b,c}. 

If xk G (ab) then cr^ G {2,4}, if xk — a then gk G {1,2,4}, and so on. For 
example, if x = (j>yr(u>iU!2^3---) and the only point on the backwards orbit of x 
which lies in (ab) U (be) U (ca) U {a,b,c} is xk G (ab) then ^ ({x}) = {cr G : 
(Tfc = ojk for all k ^ K, and ct^ G {2, 4}}. 
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Let A denote the boundary of A as a subset of M 2 , V = (ab) U (be) U (ca), and 

s = Auvu( U Mv)). 

{<7efi':M>x} 

Then 

(i) each of the sets A, V, h±(v), h 2 (v), h 3 (v), /i 4 (v), h n (v), /i 12 (v),... is dis- 
joint; 

(ii) T(S)=^(S) = H; 

(hi) T is one-to-one on h a (v) and T(h a (v)) = (v) for all a £ fi' with 

M>i; 

(iv) T is one-to-one on V and T(v) = A\{A, B, C}; 

(v) T is two-to-one on A and T(A\{A, B, C}) = T(A) = A. 

The transformation T maps A continuously onto itself. If x goes around A 
clockwise once, then T(x) goes around A anticlockwise twice. It does so in such a 
way that T{{A}) = {A}, T((Ac)) = (Ab) U {b} U (6C), T({c}) = {C} and so on. 
This information provides us with the directed graph, with labelled edges, shown 
in Figured] We denote this graph by G. It is such that that there is a bijective 
correspondence between the points of A and the set of all paths in G. A path 
in G is obtained by starting at any node and successively following edges in the 
directions of the arrows, yielding an infinite sequences of edges. The set of addresses 
of the point represented by a path in G consists of all sequences of the numbers 
{1, 2, 3, 4} which can be read off successivey from the path, with one symbol from 
each edge. For example, the only possible address for {^4} is 3 = 3333.., and the 
set of addresses of a point in (Ab) may be 3{1, 2, or 3}222..., or 33{2, 3, or 4}111.., . 
or one or more which begin 332222. 

Now suppose x £ h a (\/) for some a £ fi' with \<j\ > 1. Then by repeated 
application of (iii) above we find that the first \a\ terms in any address of x are 
precisely a, and T°'"'(x) £ V. Since T°I <T maps h a (V) one-to-one onto V, it follows 
that the set of all sets of addresses of all points in h a (v) is the same as the set of 
all sets of addresses of all points in V after a has been appended to the front of 
each of the latter addresses. So what is the set of all sets of addresses of all points 
in V? 

We have V = (ab) U (be) U (ca). Let us deal with (ab). The transformation T 
maps (ab) one-to-one onto (AB) = (Ac) U {c} U (cB). It follows that the set of all 
sets of addresses of all points in (ab), which we denote by C(ab), is determined by 
the set of sets of addresses of all points in (Ac) U {c} U (cB), which we denote by 
C((Ac) U {c} U (cB)). Specifically, 

C(ab) = {{r/a : r\ £ {2, 4}, a £ tt} : n G C((Ac) U {c} U (cB))} 

The set of sets of addresses C((Ac) U {c} U (cB)) corresponds to the set of paths in 
G which start at a node labelled (Ac), {c}, or (cB). We can similarly describe the 
address structures of (be), and (ca). We are thus able, in principle, to write down 
the set of addresses of each x £ 5; in particular, the set of all sets thus obtained 
does not depend on a, j3, or 7. 

Next we deal with S c = A\S. Since all points with multiple addresses lie in S 
and the backwards orbit of each point in S c lies in S c it follows that each point in 
S c has a unique address, and in particular <^j- 1 (S c ') = r2\^ 1 (S) does not depend 
on a, (3, or 7. 
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Figure 4. This graph can be used to compute the IFS addresses 
of all points which lie on the boundary of the triangle ▲ . 



Finally we note that = Q. Hence fljr = ft. Hence fljr = jr 1 (S) U 

(^^(S* 7 ). Since the equivence class structures of both c/)^ 1 (S) and (fi'jl^EF) do not 
depend on a, (3, or 7, it follows that equation l6.1l is true. 

So, for example, let T = ^0.5,0.5,0.5 = {>! hi, h 2 , ha, /14} and Q = ^0.65,0.3,0.4 = 
{A; gi, g 2 , 33, 34}- Then Cjf = and, by Theorem 1, the fractal transformation 
(fig o tjt : A — > A is a homeomorphism. Figure [3] illustrates the action of this 
homeomorphism. The figure on the left shows the set S, defined to be the union 
of the attractors of the two IFSs {A; hi, ha, /14} and {A; h 2 , /13, h^}. The image on 
the right shows the set S, defined to be the union of the attractors of the IFSs 
{A; gi, 33, 34} and {A; 32, 33, 34}- The two sets are related by (fig o tj?(S) = S. 
Figure [7] illustrates two other homeomorphisms associated with the family J- a ^ n . 
These examples were computed as described in section [7J 

6.4. Example 2. An example of address structures Cp and Cg such that Cjr -< Cg 
and C r ^ Cg is provided by taking T = {□; fi, f 2 , fa, U] and = {□; g x , g 2 , 33, 34} 
to be the IFSs of affine maps specified in Tables and [3] respectively. Here Del 2 
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denotes the filled square with vertices at I — (1,1), J = (1,0), K — (0,0), J = 
(0, 1). The attractor Ajr of T is represented by the fern image in Figure El The 
attractor Ag of Q is □. 

The transformations of T are such that 

(6.2) h(i) = m,h{k) = k,f 2 (i) = i,h(k) = m, 
f 3 (i) = m, f 3 {k) = I, f 4 (i) = to, f 4 (k) = j, 

where the points i,j,k,l,m G Ajr are approximately as labelled in Figure [5] Fur- 
thermore fp(Af) n f q (Ajr) — to whenever p, q G {1, 2, 3, 4} with p ^ q. It is readily 
deduced that k = <j>r(T),i = <M2), m = 0y(12) = 0^(21) = <M32) = <fcr(42), 
that 

fi r ={(re!l: 5 on (o-) ^ {21,32,42} for all n G {0,1,2,...}}, 
that Vt?r = ft and that the address structure of T is 

where 

C { r = {W} ■ a e ^,S°"(cx) ^ {12,21,32,42} for all n G {0,1,2,...}}, 
C^ 2) = {{a' 12, a' 21, a' 32, a' 42} : a 1 G O'}. 

To determine the address structure of Q, we note that □ is the union of four rect- 
angular tiles g n (0) which share portions of their boundaries. The transformations 
of Q are such that 

(6.3) 3 i (I) = M, 9l (K) = K,g 2 (I) = I,g 2 (K) = M, 
g 3 (I) = M,g 3 (K) = L,g 4 (I) = M,g 4 (K) = J, 

where the points /, J, K, L, M G Ag are approximately as labelled in Figure [5] 

Note that equations 16.31 are the same as equations 16.21 upon substitution of 
h,f2,f3,fa,hj,k,l, and to, by #1,52,53,54,^, J, K,L, and M respectively. It is 
readily deduced that K = (f>g(T),I = 4>g(2),M = (f>g(12) = <fe(21) = <fe(32) = 
0g(42), and that fig = ft. As a consequence Cjr -< Cg: if s G Cjr then cither 
s G or s G C^p ; if s G then s is a singleton and, since is a partition of 
17, there must be t G Cg such that s C t; if s G C^ then s = {a'12, o'2\, ct'32, ct'42} 
for some a' G ft' , and since M = </>g(12) = 0g(2l) = (f>g{32) = (f>g(42), it follows 
that contains a set that contains s. Hence Cjr -< Cg and, by Theorem 1, the 
fractal transformation (f>g otj: from the fern-shaped set onto □ is continuous. This 
transformation is illustrated in Figure[Hl as described at the start of section[7J Note 
however that in this case Cjr ^ Cg because there is a set in Cg which consist of a 
pair of distinct addresses, whereas all sets in C? contain either one or four distinct 
addresses. 

If, in this example, we change Q to Q specified in Table [4] then the attractor is 
still the filled square, that is Ag = □, but Equation 16.31 no longer holds and we 
can show that the fractal transformation <j)g o 7> from the fern-shaped set onto □ 
is not continuous. This lack of continuity is illustrated in Figured! as described in 
section [Jj 
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Figure 5. (i) Shows the points k, I, m and (ii) shows the points /, J, K, L, M. 



n 


a n 


bn 




d n 






1 


0.85 


-0.05 


0.125 


0.05 


0.85 


-0.039 


2 


0.06 


0.02 


0.45 


0.0 


0.165 


0.835 


3 


0.17 


0.22 


0.195 


-0.22 


0.17 


0.776 


4 


-0.17 


-0.22 


0.805 


-0.22 


0.17 


0.776 



Table 2. 



n 


a n 


bn 




d n 




1)1 


1 


0.8 


0.0 


0.0 


0.0 


0.8 


0.0 


2 


0.2 


0.0 


0.8 


0.0 


0.8 


0.2 


3 


-0.2 


0.0 


1.0 


0.0 


0.8 


0.0 


4 


0.8 


0.0 


0.0 


0.0 


-0.2 


1.0 


Table 3. 


n 


a n 


b n 




d n 




I n. 


1 


-0.8 


0.0 


0.8 


0.0 


-0.8 


0.8 


2 


-0.2 


0.0 


1.0 


0.0 


-0.2 


1.0 


3 


0.8 


0.0 


0.0 


0.0 


0.2 


0.8 


4 


0.2 


0.0 


0.8 


0.0 


0.8 


0.0 



Table 4. 



7. Pictures of tops functions 

When the underlying space is M. 2 we can use the chaos game to compute illus- 
trations of fractal transformations. 
Let two hyperbolic IFSs 

? ■= {O-Ji, ■■■Jn} and Q := {□;£/!, ...,g N } 
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and a picture function 

$)J: □-► £ 

be given, where 

□ := {(x,y) £ M 2 : < x,y < 1}. 

Let denote ^3 restricted to Ag, that is ?@g = ^\a g - Then we define a new 
picture 

by 

We say that tyjr is defined by tops plus color-stealing. 

In order to make a physical picture of tyjr and thus illustrate the tops function 
<j)g o Tjr we use a variant of the chaos game. To work at finite precision we partition 
the set DcK 2 into a finite set of small rectangles, say ten thousand of them, which 
we refer to as pixels. Each point (x, y) £ □ belongs to exactly one pixel, which we 
denote by p((x, y)). 

Start from an arbitrary pair of points (xq , y$) £ □ and (xq , j/q ) £ □. Let K be 
a large number such as ten million. For k = 1, 2, ...K let Cfc denote an element of 
{1, 2, N} chosen at random, independently of all other choices. Let 

{ x k,Vk) = f<rk( X k-liVk-l) and i x kiVk) = 9<T h ( x k-^Vk-l)- 

For each iterative step k > 100, if the color of the pixel p((x^ , y~[)) was not assigned 
at an earlier step I < k such that <7j<7/_ 1 c>7_2"-o'il > Ok (J k-\< J k-i---o\]-, then plot 
the pixel p((x^ ,Vk)) m the color ySg(x%,yj!). 

The reason this algorithm converges in practice to produce a stable physical 
picture that approximates tyg o ^ o tjt is described in Chapter 4 of (SJ. Again, it 
depends on Birkhoff's ergodic theorem. Intuitively ergodicity of the shift trans- 
formation ensures that, almost always, the sequences {(x^ , y^)} and {{xt,y?)} 
repeatedly visit all of the pixels that represent the points of Ayr and Ag respec- 
tively. Let <j( k ' = Ok crfc_i(7fc_ 2... cil. Then the point (x^,y^) is very close to 
(j>jr{(T^) when k is sufficiently large; indeed 

fo(M° {h) )M,Vk)) < ^(^(T),(^,^)). 

Similarly {xf,yf) is very close to 4>g{<J^) when k is sufficiently large. Hence, to 
a good approximation, the color of the pixel p(4>jr(a^)) is updated to become the 
color of the pixel p{4>g (a^ ')) except when > for some I < k for which 
p(M° (l) )) = P(M° (k) ))- Let 

100 < fci < k 2 < h < ... < k M < K 

denote the sequence of successive values of k at which such updates occur. Then 
{p( kl ^}fL x is an increasing sequence of addresses, each associated with a point in 
the pixel p((f>jr(a^ kl ^)). Hence, again invoking ergodicity {cr( kl '}f?_ 1 approaches the 
highest address of all points in the pixel p(^(cr^ fel ')). The address a^ kM ^ is our 
approximation to sup{rjr(a) : a £ T^(0^r(cr^ fcl ^))}. In general we expect it to 
become increasingly accurate with increasing K . According to this approximation, 
the pixel p{4>jr{a^ kl " > )) is assigned the colour of the pixel p((l>g(Tjr{a^ kM '))). Thus 
we obtain a sensible pixel-based approximation to tyg o 4>g o t?. 

In Figure [6] we illustrate two different fractal transformations from a fern-like 
fractal to a filled square, computed using this algorithm. For the picture on the left 
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Figure 6. The ferns on the left and right are both obtained by 
fractal transformations. The one on the left is continuous image of 
the central image. 



T and Q are as discussed in section T6.31 with Cjr -< Cg, Cjr Cg, so that (f>g o tjt 
is continuous. The picture tyg is represented in the center of Figure [6] It has been 
chosen to have apparently continuously varying intensity so that the continuity of 
tfigOTF is illustrated by the smooth variation of intensity in the left-hand fern image, 
which represents a close-up on *}3jr = ?Pg(4>g o tjt). To produce the picture on the 
right the IFS Q has been switched, from the one in Table [3] to the one in Table [4] 
so that <pQ ° T T is n °t continuous and ^t{<Pq ° t) 1S n0 longer smoothly varying. 

In Figure[7]we illustrate two examples, computed using the modified chaos game 
described here, in each of which the fractal transformation <pg o tjt : Ajr — » Ag is 
a homeomorphism. The homeomorphisms arc constructed using IFSs of the form 
Tu^.-y discussed in sections [2] and 16.31 In both examples Ayr = Ag = A, the 
filled triangle with vertices at A = (0,0), B = (1,0), and C = (0.5,1). Also in 
both cases, tyg : ▲ — * € corresponds to the grayscale picture of a caged bird in 
the top triangle in Figure [7] The image at bottom left shows ^g o (f>g o when 
T = ^b.525,0.525,0.525 and Q = ^0.475,0.475,0.475- In this case the corresponding 
subtriangles have the same areas at all levels with the consequence that the fractal 
transformation <j)g o is area-preserving. To produce the image at the bottom 
right we used T = ^"0.4,0.6,0.475 and Q = T . 5l 0.5, 0.5- 

8. The tops dynamical system 

In general, to determine the nature of the fractal transformation 4>g o tjf : Ayr — » 
Ag we need to know the tops code space VI jr. Here we prove that fljr is shift 
invariant. Consequently it may be described in terms of the orbits of an associated 
dynamical system Tjr : Ayr — > A jr. 

Throughout this section we assume that the transformations of the IFS T are 
one-to-one. Let Sjr : fljr — > il denote the shift transformation, defined by 

Sjr{ai(T2V3—) = cr 2 (T3a-4... 

for all oxOiV'i--- G ^f- Let 

Gjr := {(x,Tjr(xj) [ X G Ajr} 

denote the graph of the the tops function rjr. 
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Figure 7. Two examples of fractal homeomorphisms applied to 
the picture at the top. The tranformations from the top image to 
the one at bottom left is area-preserving. 



Lemma 3. Let (x,a) £ G r . Then (f- i 1 (x),S r (a)) £ G r . 

Proof, (x, a) £ Gjr implies x £ Ayr, a £ ftp and Tjr{x) — a. In particular, 

<j)jr[a) = x for any z £ X, 

lim f ax o / ff2 o ...fa k {z) = x. 

k — >oo 

Using the continuity and invertibility of f ai it follows that 
lim f a2 o f a3 o ...fa k {z) = f~ X {x). 

This says that <t>jr{Sjr{a)) = f-?{x) which tells us that S r {a) £ 4^ X {{f^{x)}). 

Now suppose that there is w £ 0^ 1 ({/<7 1 1 (a;)}) with u> > Sp{o). Then (j>r(u>) = 
fa X { x ) which implies f ai (ipyr[u>)) = 4>f{o~\lo) = x. Let a = o\lo. Then a > a and 
4>r(a) = x which contradicts the assertion that a is the largest element of il such 
that (j>j:{o-) = x. Hence Sp{a) £ VLjr and Tjr{f~^{x)) = Sp{a). □ 

Lemma 4. Let (x,a) £ Gjr. Then (fi(x),la) £ Gjr. 

Proof, (x, a) £ Gj: implies tjt(x) = a. Hence x — 4>f((t) and so fi(x) = (j)jr{\o~). 

Now suppose that (fi(x), la) ^ Gjr. Then there is u> > la such that 4>f{uS) = 
fi(x). But then lo = la where a > a and (j)^(la) = fi(x). This implies 4>r(a) = x 
with a > a which implies tj^(x) > a which is a contradiction. Hence (fi(x), la) £ 

Gjr. □ 

It follows from Lemmas [3] and [4] that the mapping T? : Gf — > Gf specified by 

%{x,a) = (f-\x),Sr(a)) for all (x,a) £ G T 

is well-defined and onto. 

In particular, the projection of Tjr on fijr yields the symbolic dynamical system 
Sp : Qjr — > because from Lemma |4] we have 
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The projection of Tjr : Gjr — > Gjr onto Ajr yields what we call the tops dynamical 
system 

Tjr \ Ajr — ► Ajr 

where 

f fx\x) if xeD 1 :=fi{Ar), 

fc\x) if x£D 2 :=f2(A F )\f x (A F ), 



.1) T r (x) = 



N-l 



fx\x) if x£D N := ftr(4F)\ U /t.(Af), 

n=l 

for all x S Ajr. Lemma @] implies 

Tjr(Ajr) = A jr. 

Theorem 2. TTie tops dynamical systems Tj= : A? —> Aj= is related to the symbolic 
dynamical system Sjr : f2jr — > £7jr 6 2/ i/ie tops function Tjr : Ayr — » J7jr 7 according to 

Sjr = Tjr oTjr . 

If ftp C fig then 

(4>g o tjt o 2»(ar) = (T g o g o r^)(a;) /or a// a; £ Ajr. 

J/ Cjf = Cg iften toe tops dynamical systems Tjr : Ajr — > Ajr and Tg : Ag — > Ag ore 
topologically conjugate. 

Proof. Let <£>jr = r^ 1 be as discussed at the end of section 2J Then we claim that 

Tjr O $jr = $jr o ,Sjr. 

Since <5>jr maps iljr onto itself and <I>jr maps iljr onto Ajr it follows that the mapping 
<I>jr o Sjr takes fl jr onto Ajr. (Similarly, Tjr o $jr maps fl jr onto A jr.) 
Let a = <7io~20"3-" G flp- Then iSjr(<7) = CT2C3... £ fijr and 

(4>jr o5jr)(a) = lim (f a2 o / CT3 o ... o / CT J(z). 

A; — >oo 

On the other hand 

$jr(o-) = lim (/ CT1 o / CT2 o ... o /„„)(«) 

ft — »00 

= UA lim (/ CT2 o / CT3 o ... o f ak )(z)) = ($^(0-203...)), 

& — >oo 

belongs to Ajr and lies in the range of / CT1 and so must belong to D ai as defined in 
Eauation l8.il Hence 

(Tjr o $jr)(cr) = $jr(cr 2 cr 3 ---) = ° 5jr)(cr) 

for all a € fljr. 

We now apply Tjr to both sides of this last equation to complete the proof of the 
first assertion in the theorem. 

Now assume that fijr C fig. Then, since 

Sjr-(o-) = Sg(a) 

for all a £ fijr, it follows from the first part of the theorem that 

(Tjr o Tjr o $jr) (a) = (rg oTg o $g)(cr) 
for all a £ 51 jr. It follows that 

(Tjr o Tjr o $jr o Tjr)(x) = (rg o Tg o $g o Tjr){x) 
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for all x € Ayr. But <&jf ° tjf = and 

($g O T^)(X) = (0g O Tyr)(x) 

for all x e Ajf. Hence 

(TjF O Tyr){x) = (Tg oTg Olfig O Tyr)(x) 

for all x S Ajf- Applying </>c; to both sides we obtain 

(4>g O Tyr O Tyr) (x) = ((fig O Tg O Tg O (j)g O Tyr) ( X ) 

for all x G Ajf. But (j)gorg = iA g - This completes the proof of the second assertion 
in the theorem. 

Finally, let us suppose that Cyr = Cg. Then Theorem 1 implies that <f>g o Tyr is a 
homeomorphism from Ayr onto Ag. Also Cjf = implies Qyr = Qg which implies, 
via the previously proven part of this theorem, 

Tyr(x) = {(j>g O Tjr)^ 1 O Tg O (f>g O Tjr)(x) 

for all x e Ayr. □ 

If the domains {D n : n = 1,2, .., N} are known then it is easy to compute the 
tops function. Just follow the orbit of x under the tops dynamical system and keep 
track of the sequence of indices ci 02 03 04... visited by the orbit. 

In the special case where the IFS is totally disconnected and the /„s are one- 
to-one then T : A — > A is defined by T(x) = fn ( x ) where n is the unique index 
such that x 6 f n (A). This dynamical system has been considered elsewhere, for 
example in [l] and [14] . In this case 4>f • fl — > Ayr is a homeomorphism, Tyr = <f>^ , 
and Tyr it is conjugate to the shift transformation according toTjr = <pyr o Syr o 4>~yr X ■ 

Theorem [5] says in particular that Tyr : Ayr — » Ayr is a factor of Syr : Vlyr — > Q,yr, 
and as defined for example in [T3] p. 68, because QyroSyr = Tyro&yr where $>yr = r^ 1 
is continuous; this tells us that the topological entropy of Tyr is less than or equal to 
the topological entropy of Syr, [13] Proposition 3.1.6, p. 111. If Tyr is continuous then 
Theorem [2] says that the two dynamical systems Tyr : Ayr — > Ayr and Syr : ilyr — > VLyr 
are topologically conjugate, see [13] p. 60, and it follows that the two systems must 
have the same topological entropy. 

This suggests that we may compare the complexity of some subsets of R 2 by as- 
signing to them the topological entropy of a corresponding shift dynamical system. 
Let M denote the set of all attractors of hyperbolic IFSs in M 2 , whose trans- 
formations are all affine and invertible, such that the associated tops function is 
continuous. Then we can define the topological entropy of each Ayr to be the infi- 
mum of the entropies of the set of corresponding shift dynamical systems. In this 
way we arrive at a geometry-based definition of the topological entropy of some 
subsets of M. 2 . Is it useful? 
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